Intermediate Coupling Theory of Electronic Ferroelectricity 
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We calculate the quantum phase diagram of an extended Falicov-Kimball model for one and two-dimensional 
systems in the intermediate coupling regime. Even though some features of the phase diagram are obtained 
analytically, the main results are calculated with a constrained path Monte Carlo technique. We find that this 
regime is dominated by a Bose-Einstein condensation of excitons with a built-in electric polarization. The 
inclusion of a finite hybridization between the bands removes the condensate but reinforces the ferroelectricity. 
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In a recent paper jjl, one of us demonstrated that a novel 
ferroelectric state, originally proposed by Portengen ef a/ 
is present in the strong-coupling and mixed-valence regime of 
an extended Falicov-Kimball model (FKM) |S]. In contrast 
to the traditional ferroelectric (FE) transitions, which are a 
subgroup of the structural phase transitions, the new FE state 
is caused by a purely electronic mechanism of an essential 
quantum nature. Its electronic origin and the strong interplay 
between the orbital and the spin flavors of the electronic de- 
grees of freedom make this novel state suitable for new tech- 
nological applications. For example, it opens the possibilities 
of faster switching ferroelectrics and controlling their optical 
properties with magnetic fields. Accordingly, it is important to 
define more clearly the conditions for the experimental obser- 
vation of this novel state, especially in light of several recent 
experiments that might have seen it 0,01 ■ 

Our extended FKM model consists of two dispersive bands 
plus an inter-band local Coulomb interaction. In absence of 
explicit hybridization between the bands, the half-filled model 
exhibits a FE state that is a Bose-Einstein (BE) condensate of 
excitons. The electric polarization is the result of a sponta- 
neous phase-coherent hybridization between the two local or- 
bitals with opposite parity under inversion. The FE moment 
is proportional to the real part of the U{\) order parameter 
of the condensate. These conclusions were based on a strong 
coupling treatment of the model. An accompanying perturba- 
tive analysis suggested that if hybridization between the bands 
were included, the BE condensate (BEC) would be replaced 
by a distinguishing Ising-like FE state Because hybridiza- 
tion is natural and the experimental consequences are signifi- 
cant, resolving this point is important and timely. 

To resolve this point, we will now step way beyond the orig- 
inal work by moving away from the strong-coupling limit and 
from the perturbative treatment of the hybridization and fo- 
cus on the intermediate coupling regime. We find that ferro- 
electricity remains robust. In particular, it possesses the same 
phases as the strong and weak couplings regimes, even though 
the nature of the condenstate changes continuously from BCS- 
like in the weak coupling limit to a BEC of hard-core bosons 
in the strong couping regime |6|||. We will in fact report the 



complete one (ID) and two dimensional (2D) quantum phase 
diagrams of the extended FKM in the intermediate coupling 
regime. The results are obtained from a scaling analysis of 
the numerical data produced with the Constrained Path Monte 
Carlo (CPMC) technique |3|. When we include a finite hy- 
bridization between the bands, we find as predicted by the per- 
turbation analysis a very large FE (or antiferroelectric) region 
for the mixed valence regime which is no longer a condensate 
of excitons but a state with an Ising-like symmetry that has 
a built-in electric polarization. Importantly, the intermediate 
coupling phase diagram indicates more accessible regions of 
the ferroelectric phase than in the strong and weak coupling 
regimes. 

There are two recent experiments to which our results might 
apply. Moro et al observed a FE state in free Nb clusters f^l- 
The experiment shows that cold clusters may attain an anoma- 
lous component with very large electric dipole moments. It is 
clear that in pure Nb the FE moment cannot be produced by 
a lattice distortion. Therefore, the ferroelectricity must have 
an electronic origin. In addition, ab initio calculations fllll 
in pure Nb indicate that the valence electron has the same 
amount of s and p character This observation makes the 
present theory a candidate to explain the anomalous ferroelec- 
tricity of the free Nb clusters. 

Another relevant experiment is the phase coherent state ob- 
served by Spielman et al | 5]. The two bands are provided by 
two parallel 2D electron systems interacting through an inter- 
layer Coulomb interaction. There is also an applied magnetic 
field whose value is such that the Landau level filling factor of 
each layer is 1/2. From the results obtained in Ref. Q]] and 
the theory that we include below, we propose an experiment 
to observe electronic ferroelectriciy in these bilayer systems. 

Our extended FKM for spin-less fermions on a D- 
dimensional hypercubic lattice is: 

H = ^6.<+ *-'(/iUi+e„.'+//+e,.'/i.) 

+ C/'^''^K-l/2)(nf-l/2), (1) 

i 

where = {a, b} is the orbital index, e^^ is a vector in the i] 
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direction (ij — {xi, X2, xd}), and = hufiv '■^^ occu- 
pation number of each orbital. Note that tab = —ha because 
the two orbitals have opposite parity under spatial inversion. 

This Hamiltonian can be rewritten as an asymmetric Hub- 
bard model if the orbital flavor is represented by a pseu- 
dospin variable = f^^, = /;^, Cj^ = f^^^ and 

C;| = The pseudospin components are given by: r" = 

5 J2uy fL^""' f\v'^ with a {a;, y, z}. The new expression 
for H is an asymmetric Hubbard model with a Zeeman term: 

i,(T (iJ),o-,cr' 

+ ^nitnu+S.^r/, (2) 

i i 

where = ^(ca + ttT = ^qq' = hb, ^Ti *Qb' 
^it ~ ha and Bz = ea — If there is no hybridization be- 
tween the bands, there is a remaining U{\) symmetry gener- 
ated by the total magnetization along the z axis ]\F- = t(' . 
In the original language, this U{1) symmetry corresponds to 
the conservation of the difference between the total number 
of particles in each band (the bands cannot interchange parti- 
cles). 

We will consider from now on only the half-filled case, i.e., 
one particle per site. The ground state of H exhibits two pos- 
sible magnetic orderings in the strong coupling limit: Ising- 
like antiferromagnetism and xy-like ferromagnetism (antifer- 
romagnetism) for 7 = — ij,i/^TT > (7 < 0) (Ol- In terms 
of the original variables, the Ising-like phase corresponds to 
staggered orbital ordering (SOO) and for D > 2, the xy-like 
ferromagnetism corresponds to a BEC of excitons with a built- 
in electric polarization Q]: 

P=i]E(/ia/i6 + /ia//j = |^M-, (3) 

where fi is the inter-band dipole matrix element and 51 is the 
volume of the system. The electrical polarization is then pro- 
portional the real part of the U{1) order parameter of the con- 
densate. 

An addtional particle-hole transformation on one of the spin 
flavors maps Eq. |2] into a negative U asymmetric Hubbard 
model with the Zeeman term replaced by a chemical potential 
01 ■ Under this transformation, the xy-\ike ferromagnetism is 
mapped into a superconducting state. This establishes a for- 
mal connection between the BEC of excitons and a supercon- 
ductor. An eventual condensate of electron-hole pairs in the 
weak coupling regime is then formally analogous to a BCS- 
like superconductor 

One dimensional case. From now on, we will set 7 > 
and t|| as the unit of energy. We start by on considering 
the simplest case of zero hybridization between the bands: 
t|l = = 0. In Fig.^ we present the ID quantum phase 
diagram of as a function of 7 and Bz ■ Most of the phase 
diagram was calculated with the unbiased CPMC technique 
0]. The diagram shows the same three phases which were 
obtained in the strong coupling limit 12]. For Bz = 0, the 



ground state is an Ising-like antiferromagnet (SOO). At a criti- 
cal value of I = B'f^, there is a transition from the orbitally 
ordered gaped phase to a critical (gapless) phase with dom- 
inating pseudospin-pseudospin correlations in the xy plane. 
This second phase is also characterized by a non-zero uniform 
magnetization along the z axis. When \Bz\ reaches a sec- 
ond critical value, _B^^, Af^ saturates and a new gap opens. 
In terms of our original language, this is the transition from 
a mixed valence to non-mixed valence regime in which one 
band is full and the other one is empty. This is also a transi- 
tion from an excitonic to a band insulator 
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FIG. 1: (color online) One dimensional quantum phase diagram of 
the extended FKM U"''' = l{tti = = 0) 01. The circles are the 
calculated points for BJ^ and the full line is an exponential fit. The 
error bars are smaller than the symbol sizes. The squares are numer- 
ical results for B?^ and the full line is the analytical calculation. 

The value of is exactly obtained by solving the one 
particle problem of flipping one pseudospin in a fully sat- 
urated system. To show that this value is not altered by a 
discontinuous change of the magnetization, we also include 
the numerical calculation of B'f based on the computation of 
£'m (A/^, B^ ~ 0), minimum energy for the subspace with a 
given magnetization AI^. Note that this can be done because 
AI^ is a good quantum number of H{t^i = i|f = 0). is 
the critical value of B^ which is required to change the total 
magnetization Af^ of the ground state, i.e., it is the minimum 
value of [E,niM^,B^ = 0) - Egs{B^ = 0)]/A/^ where 
Egs is the ground state energy. To calculate in the ther- 
modynamic limit, we scaled this minimum value for different 
chain sizes up to Ns = 100 {Ng is the number of sites). A 
typical scaling is shown in Fig.|2^ for two different values of 
7. The finite extrapolated value of i?^^ implies a finite gap A/ 
for the Ising-like (SOO) phase. The case 7=1 corresponds 
to the SU{2) critical point for which the pseudospin correla- 
tions along the three different axis are the same. For the strong 
coupling limit Q], it is known that the gap of the Ising-like 
(SOO) phase closes exponentially (Kosterltiz-Thouless tran- 
sition) when 7 gets close to one (see Fig.Q. 

The different phases were identified by computing the 
pseudospin-pseudospin correlation function: 

^"(q) = E(^"^" - (Tr)(^"))e^'>-("-'-'-^) (4) 
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In Fig.|2j5, we show the scaling of S^' {0) and S'^(7r) calculated 
with the CPMC for = and 7 = 0.1. As expected, S'^in) 
has a finite value in the thermodynamic limit indicating that 
the ground state has long-range staggered orbital ordering. 
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FIG. 2: (color online) U=l, f|j = — 0. a) Scaling of calcu- 
lated with the CPMC method in chains of length Ns. b) Scaling of 
S^(0) and 5"'(7r) for — 0. c) Scaling of the energy gap A^y for 
/Ns = 1/4. The scaling function is a second order polynomial 
expression in 1 /Ns . d) Pseudospin-pseudospin correlation function 
for M^Ns = 1/4. 

In Fig. |2|;, we show the scaling of the energy gap of the 
xy-like phase, A^^j,, for /Ng = 1/4. A^y scales to zero in 
the thermodynamic limit as expected from the critical nature 
of this phase. Fig.lJJl shows that the pseudospin-pseudospin 
correlation function (see Eq. |4} in the x (or y) direction is 
much larger than the corresponding correlation along the z 
axis. This shows the xy-like nature of this phase. 

Two dimensional case. The main qualitative difference be- 
tween the ID (Fig[0 and the 2D (Fig. |3} quantum phase dia- 
grams of H{t-^i = = 0) is that the excitonic condensate is 
critical in one dimension, as required by the Mermin- Wagner 
theorem fl2il . and has long range order in two dimensions. In 
the pseudospin language, this means that the 2D system has 
a finite in-plane magnetization (see Fig.|3. In addition, 
Eq. Q shows that this state has built-in electric polarization 
proportional to the x component of or the real part of the 
complex order parameter of the excitonic condensate Q). 

As in the ID case, B'fij) was computed exactly, and the 
transition from the SOO to the ferroelectric BEC of excitons is 
simultaneously a first order valence transition because the rel- 
ative population between the two bands, M^, changes discon- 
tinuously. Note that in the intermediate coupling regime, the 
extension of the ferroelectric phase for 7 not very close to zero 
is of the order of the bandwidth while in the strong coupling 
regime the size of this region is proportional to taatbb/U""^. 

We conjecture that the three dimensional (3D) quantum 
phase diagram of H should be similar to the 2D one since 
the quantum fluctuations are smaller in the first case. More- 
over, we expect the quantum phase diagram to have the same 
topology for the weak and intermediate coupling regimes in 
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FIG. 3: (color online) Two dimensional phase diagram of the ex- 
tended FKM for U = 2 (tn = tif = 0).The circles are the calcu- 
lated points for B'^^ and the full line is a guide to the eye. The error 
bars are smaller than the symbol sizes. The squares are numerical 
results for and the full line is the analytical calculation. 



any dimension D > 2. The main expected difference between 
D=2 and D=3 occurs at finite temperatures due to the Mermin- 
Wagner theorem |12J. The transition temperature associated 
to the BE condensation is finite only for the 3D case. The BEC 
of electron-hole pairs undergoes a Kosterlitz-Thouless phase 
transition in a 2D system. 

From Eq. Q, we see that in any dimension the observed 
electronic hybridization is the consequence of a spontaneous 
phase-coherent hybridization between the two bands. The 
phase must have a real component in order to produce a fi- 
nite electric polarization. We can ask what is the physical 
consequence of a pure imaginary phase cj) = ±7r/2 (tan0 = 
/Ajy), i.e., a magnetization pointing in the y direction. In 
this case, instead of ferroelectricity, the phase-coherent state 
induces an ordering of atomic cuiTents |'13>]. Since the two 
atomic orbitals have opposite parity, these currents do not 
produce a net local magnetic dipole but they can have a net 
quadrupolar magnetic moment. Therefore, the y direction in 
the pseudospin space corresponds to ferroquadrupolar or anti- 
ferroquadrupolar magnetic ordering. 

As predicted in |1], we will now show that the inclusion 
of an explicit hybridization between the bands breaks explic- 
itly the C/(l) symmetry which is spontaneously broken in the 
condensate and removes the degeneracy associated with all 
the possible orientations of M^. Along with this the natu- 
ral question to answer is: what is the particular orientation 
favored by a non-zero hybridization. 

Inclusion of Hybridization. In addition to the trivial U{1) 
symmetry associated with the conservation of the total num- 
ber of particles, the remaining symmetries of H for non-zero 
hybridization are the spatial (/) and the temporal (T) inver- 
sions. Since the two orbitals have opposite parities, the signs 
of and are interchanged under the application of /. 
For the pseudospin components we have: /rf / = —Tf , 
It^I = — ty' and It[I = t^. The temporal inversion (com- 
plex conjugation) leaves rf^ and t? invariant and changes the 
sign of T^. The ferroelectric state occurs when / is sponta- 



3 



LA-UR-OO-XXXX 



February 6, 2008 



submitted to Physical Review Letters 



neously broken. The ordering of atomic currents occurs if 
both / and T are spontaneously broken. Note that the product 
IT is still a symmetry of this second state. 
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FIG. 4: (color online) Scaling of S^XO) and S'^(O) for a 2D system: 
[7 = 4, t^i = = 0.2 and = 1 (full lines), (7 = 2, 

hi = -iiT = and W /Ns = 1/4 (dashed line). The scaling 
function is a second order polynomial expression in 1/Ns. 

In Fig. m we plot S'^(O) and 5^(0) as a function of the 
size of a two-dimensional system. The largest system has 
Ns = 100 sites. The scaling clearly shows that the ferroelec- 
tric state {Mx ^ 0) is stabilized in the thermodynamic limit. 
This agrees with a perturbative analysis which is only valid in 
the strong coupling limit Q] . 

The BEC of excitons is thus replaced by an Ising-like fer- 
roelectric state when a finite hybridization is included. This 
result has important physical consequences since most real 
systems have an explicit hybridization between two bands. 
This hybridization is the reason why a BEC of excitons is 
so hard to observe. The only exception is the bilayer system 
that we mentioned in the introductionj^l ■ According to the 
theories presented here and in Ref. [l|], this bilayer system 
must exhibit the phenomenon of electronic ferroelectricity if 
the two layers are made of bands with opposite parity (like 
"s" and "p") instead of being identical. Then, the d\P\/dE 
response as a function of the uniform electric field E should 
show the same resonant behavior observed in the tunneling 



conductance as a function of the interlayer voltage Ql • 

In summary, we derived the quantum ID and 2D phase 
diagrams of the extended FKM in the intermediate coupling 
regime. In the absence of hybridization, we found that the in- 
sulating phase obtained at half filling has a transition from a 
non-mixed valence to a mixed valence regime as a function of 
the energy difference between the centers of both bands. At 
this transition the system changes from a band to an excitonic 
insulator which is a BEC of electron-hole pairs. Addition- 
ally, we showed that this condensate contains two different 
orderings which are degenerate: ferroelectric order and chi- 
ral order. The chiral state is a spontaneous ordering of local 
atomic currents. Since the two bands have opposite parity un- 
der spatial inversion, these currents can only produce a non- 
zero quadrupolar magnetic moment. For bands with the same 
parity, the chiral state can produce orbital magnetism and the 
ferroelectric state is replaced by quadrupolar electric ordering. 
If the distance between both bands is further decreased, there 
is a first order valence transition from the BEC of excitons to 
a state with staggered orbital ordering. 

Addressing an issue most important for physical systems, 
we also considered the effect of a finite hybridization on these 
ordered degenerate phases. We showed that the U{1) degen- 
eracy associated with the BEC is lifted by the hybridization, 
and the degenerate ground state is replaced by an Ising-like 
EE state (broken Z2 symmetry). This result is valid even be- 
yond the perturbative regime jtabl << I ^aa 1, 1^66 1- It is crucial 
for applications of the present theory to systems like the free 
Nb clusters fj] which are most likely in the intermediate cou- 
pling regime and have a considerable hybridization between 
the s and the p valence bands. 
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